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ABSTRACT 

This paper describes a variable-gain, optimal, discrete, 
output feedback design approach that is applied to a 
nonlinear flight regime. The flight regime covers a wide 
angle -of-attack range that includes stall and post stall. The 
paper includes brief descriptions of the variable-gain 
formulation, the discrete-control structure and flight 
equations used to apply the design approach, and the high 
performance airplane model used in the application. Both 
linear and nonlinear analysis are shown for a longitudinal 

four-model design case with angles of attack of 5°, 15°, 
35°. and 60°. Linear and nonlinear simulations arc 

compared for a single-point longitudinal design at 60° 
angle of attack. Nonlinear simulations for the four-model, 
multi-mode, variable-gain design include a longitudinal 
pitch-up and pitch-down maneuver and high anglc-of- 
attack regulation during a lateral maneuver. 

INTRODUCTION 

High angle-of-attack (a) flight is a desirable capability 
because of potentially large payoffs for combat aircraft. 
These airplanes will have to operate over a wide flight 
regime that includes stall and post stall. Post-stall 
technology (PST) fighters must therefore be capable of 
controlled flight beyond the maximum aerodynamic lift 
angle of attack. One goal is to make airplanes roll faster at 

high a, which drives up the pitch and yaw control power 
requirements. The problem is that the control power of 
conventional airplanes decreases as airspeed decreases. 
Research is under way on methods to improve 

controllability at high a where control power is low and 
departures most likely. One solution to the control power 
problem is the use of thrust vectoring controls. The use of 
thrust vectoring provides additional control power for 
stability augmentation as well as for maneuver 
enhancement in the stall and post-stall flight regimes. 

The requirement for airplanes to operate with agility 
over a highly nonlinear flight regime puts additional burden 
on control law designers. A linear single-point design with 
constant feedback gains is probably inadequate to maintain 

good control properties over the complete a range. 
Typically, control designers develop satisfactory control 
laws at several points over the flight regime and then use 
an interpolation technique, often a least squares fit, to 
obtain final control gains. One potential problem with this 
approach is that the controller may lose performance 
characteristics and even stability, in high-order, highly- 
sensitive plants. 

An objective of this research is to extend the operating 
range of the control law over the flight regime while 
continuing to use established linear control design and 
analysis techniques. In other words, the objective is to 
design a nonlinear control law using linear theory. The 
control methodology is an extension of previous 


developments for a direct-digital feedback design 1 . This 
paper describes the application of a variable-gain optimal 
output feedback control design methodology 2 ' 4 in which 
the feedback gains are calculated and scheduled as a 

function of a. In general, several design parameters, 
cither linear or nonlinear, may be used. Desired control 
characteristics are specified prior to designing the 
controller, negating the need for interpolation or gain- 
fltling techniques. Each of the operating points considered 
in the design contribute to the minimization of a global cost 
and arc guaranteed to be stabilized. 

The variable-gain formulation is discussed in the first 
section of this paper. Following this section is a 
description of a discrete control law structure that has been 
used successfully in other direct-digital design 
applications 5-8. The third section includes a brief 
description of the airplane model used for this application. 
A section on the design example includes a description of 
four longitudinal linear design models and their associated 
linear analysis. In addition, linear time analysis for a 
single-point design is shown for comparison with 
nonlinear simulation results. The section on nonlinear 
simulation includes some practical aspects that must be 
considered, as well as single-point results and variable- 
gain scheduling results. Two different maneuvers arc 
shown for the variable-gain scheduling example. 

VARIABLE-GAIN FORMULATION 

Variable-gain output feedback is an outgrowth of a 
stochastic, optimal, discrete, output feedback design 
approach for single-point control designs 1 . Typically, 
single-point designs are combined by using a curve fitting 
approach, such as a linear approximation or a higher-order 
least-squares fit. In contrast, the variable-gain 
approach 2 ' 4 solves a control design problem by 
minimizing a performance index for an ensemble of 
systems. Equations describing the stochastic discrete 
output feedback system arc functions of the argument , p, 
which in the general case represents various operating 
point parameters. The dynamic system, measurement and 
control law equations arc 

x(p,k+l) = d>(p) x(p,k) + F(p) u(p,k) +w(p,k) (1) 

y(p,k) = C(p) x(p,k) + T|(p,k) (2) 

u(p.k) = -K(p) y(p,k) (3) 

where <t>(p), f'(p) and C(p) represent the system matrices, 
x(p,k), y(p,k), and u(p,k) are the state, measurement, and 
control vectors, w (p,k) and q(p,k) represent state process 
noise and sensor noise vectors, K(p) is the feedback gain 
matrix and k is the sampling step. It is assumed that 

w(p,k) and rj(p,k) arc uncorrclated, zero-mean, white- 
noise vectors as 


E[\v(k,p)] = 0, E[Tl(k,p)] = 0 ( 4) 

E[w(k,p) wT(j, p)J = W(p)8 k _. (5) 

Et(Tl(k,p)Tl T 0,p)] = V(p)8 kj (6) 

where W(p) and V(p) represent process-noise and sensor- 
noise covariance matrices, 8 fc . is the Kroneckcr delta and 
E is the expectation operator. 


within each control channel. Only measurable states are 
being fed back and selected outputs are being integrated to 
give type-1 control properties. 


iiu, nr equauons ana weightings for the local cost 
f“" c 4° n formulated in the continuous domain and then 
transformed to the discrete domain. The discrete plant 

equations 5 as" 0) augmcnled b 7 additional dynamic 


The feedback gain matrix has a linear functional 
relationship with the operating point parameters^ 


K(P) = K 0 + £ p. K i 
i=l 


(7) 


It is not feasible to compute the gain K(p) for every 
possible variable (a, (3, etc.) over an entire night regime, 
so a selected set of p. is used for the variable feedback-gain 

schedule. This selected set of parameters will be referred 
to as the variable-gain scheduled parameters in the 
remainder of this paper. Each of the p | can be either linear 

or nonlinear, even though (7) shows a linear functional 
relationship between feedback gain and the gain-scheduled 
parameters. In the example design problem described in 
this paper, one gain-scheduled parameter that is a linear 
function of angle of attack (a), is used. 

A quadratic performance index for the local discrete 

cost J(p, K(p)), of each plant described in (1) and (2) has 
the form 

1 N 

J(p,K(p)) = lim lE[X(p,k + l) T Q(p)X(p,k+l) + 

N— m» v 7 k=0 

U(p,k) T R(p)U(p,k)] (g) 

* 0 ’ R (P) * 0 and X(p) and U(p) (to be 
described later) are augmented stale and control vectors 
respectively. A global cost is expressed by 


M 


J(K) = £fjj (Pj.KfPj)), f>0 


u(p,k+l) = u(p,k) + (AT) v(p,k) (10 ) 

z(p,k+ 1) = z(p,k) + (AT)H z x(p,k)+ (AT)D z u(p,k) (1 1) 

where H z and D z are transmission matrices for the 
integrator states, z(p,k), v(p,k) is the control rate command 
and is used in equation (8) as U(p,k), and AT is the 
discrete sampling period. The PIF state vector X(d k) 
and output vector, Y(p,k), arc defined as P 


X(P,k) = fx(p,k) T , u(p,k) T , z(p,k) T ]T 
Y (Pik) = fy(p,k) T , u(p,k) T , z(p,k)T]T 


( 12 ) 

03) 


(9) 


equation (^ X(P,k) ' S then USed in the local cosl 

The command generator uses an output model- 
following approach to achieve specific objectives, such as 
changes in pitch rate or angle of attack. Feedforward 
gains, connecting the command model to the control 
arc . ca,c “* atcd open loop, independent of the 
feedback gains. Thc objective of CGT is to cause selected 
outputs, y, of the airplane to track the output, y , of a 

linear command model using feedforward control dieory. 

more complete description is presented in reference 10 
and chapter 4. 10 of reference 9. 

All equations are implemented in incremental form 
using total measurable quantities. This implementation has 
worked very well with nonlinear dynamics in other 
applications 5 - 7 . A significant advantage of this 
implementation is that it performs automatic trim The 
airplane goes to a new equilibrium state as a function of 
slowly varying commands and, therefore, trim values are 
not needed m flight. 


where each of the local costs arc summed and weighted 
j - t0 assign relative priorities to the individual plant 
models. 

The main point of this variable-gain formulation is to 
find a variable feedback gain matrix, K(p), which 
minimizes the global cost (9) subject to the constraints 
defined in (1) -(7). 

CONT ROL LAW structure 

A PIF/CGT (proportional-intcgral-filter with command 
generator tracking) discrete controller structure is used for 
the control lawS.9. pip i s a direct digital integrated 
lormulation using linear dynamics for design The 
controller inputs represent rate commands and the time step 
from rate to control position commands accommodate 
necessary computation delays. Proportional feedback of 
control position commands allows additional filtering 


The main incremental flight control equations^ are 

V u m,k+I ' "m, k' (14) 

Vr \ 4 (A1 > k G5) 

where y z is the vector of integrated outputs, y is the 
output vector from the command model, u is the input 
vector to the command model, K y is the feedback gain 
matrix associated with the plant outputs, and T , T-, and 
T f re P r esent precomputed matrices. The control" position 
command to the plant servo actuators is the vector u in 
equation (15). 


? 



AIRPLANE MODEL 

The model used for this work has the characteristics of 
a high performance, high angle-of-attack (a) airplane that 
can undertake fighter or attack missions. For this study, 
the airplane has a gross weight of approximately 
33,000 pounds, a wing span including tip missiles of 40 
feet, and the length is 56 feet. Controls include two 
afterburner engines plus the following control surfaces (1) 
horizontal stabilators, (2) full-span leading-edge flaps, (3) 
trailing-cdge flaps, (4) ailerons, and (5) twin vertical 
stabilizers. In addition, pitch and yaw thrust vectoring 
controls have been added for both longitudinal and lateral- 
directional maneuvering. 

A nonlinear six dcgrcc-of-frecdom, rigid-body 
dynamic model is used for batch simulations. The 
aerodynamic tables, generated from a wind tunnel derived 

data base, have an a range from -10° to 90° and a sideslip 

range of ±20°. Most coefficients are determined by using 
table lookup with linear interpolation. Flexibility effects 
are incorporated using flex/rigid ratios. During simulation, 
the aerodynamic control surfaces are driven by 
approximate transfer functions with constant rate and 
position limiting. Thrust vectoring actuator dynamics 
consist of a simple first-order low-pass filter for all four 
vanes, two for each engine, that provide pitch and yaw 
moments. Afterburner dynamics arc included in the engine 
model, and engine thrust levels based upon airflow rate, 
which is determined by table lookup. 

Table I shows the short period eigenvalues and trim 
Mach numbers for an a range from 5° to 60° at 
15,000 feet. The airplane is trimmed straight and level up 

to a equal to 35 °; as a becomes larger, the flight path 
angle becomes negative. The reason is insul 1 icicnt thrust 
is available at a greater than 35°. The data in Table 1 show 
that the short period open-loop eigenvalue becomes very 
lightly damped when a exceeds 40° (post-stall flight 
regime) with the airplane being slightly unstable for a at 
50° and 55°. The phugoid (not shown) becomes more 
heavily damped as a increases. 

The example application described in this paper 
includes the horizontal stabilator and thrust vectoring 
controls. The linear design assumes fourth-order 
dynamics for the stabilator (poles at -14.9 ± j33 and -62 ± 
j85) and first-order dynamics for the thrust vectoring 
control (pole at -30). 

RESIGN EXAMPLE 
I. inear Design Modd 

Four flight conditions (a = 5°, 15°, 35°, and 60°) arc 
used in the variable-gain output feedback application. The 
control law at these four flight conditions include two 
command modes with a transition region between these 

modes. A pitch rate, q-command, mode is used at low a 

(5°), whereas an a-command mode is used when a is 35° 
or higher. A transition region occurs between these two 
modes. The a-command is phased in between 5° and 15° 
of a and the q-command is phased out between 15 J and 
35° of a. This transition range may not be the optimum 


choice, but the example docs demonstrate a practical 
capability to incorporate multi-mode design. 

The airplane model used in the design example is 
discussed in the previous section. The plant has fourth- 
order longitudinal dynamics for the short period and 
phugoid modes, two controls with fifth-order actuator 

dynamics for the horizontal stabilator ,5 S , and symmetric 

thrust vectoring ,8 v , as well as three outputs for q, a, and 
n (normal acceleration). The n z measurement is 

transmitted through a low-pass filter of 50 rad/sec, in order 
to make the measurement a function of one state as 
opposed to several stales and controls. This additional 
filter gives a total of 10 states for the design plant. 

The PIF controller used in this example has six states. 
Two of the states are created by feeding back the control 
position commands, which generates additional low-pass 
filtering. Two additional slates are created by adding 

integrators to either q or a (depending on the mode) and to 
8 The q or a integrator provides zero steady-state error 

for the commanded variable, while the 8 y integrator 

allows the thrust vectoring control to be maintained near 
zero-deflection or some other reference input. The final 
two controller states are created by second-order dynamic 
compensation to provide additional regulation of the short 
period mode. This compensation is designed by heavily 
weighting the error between the two dynamic compensator 
state outputs and the open-loop short period modern. 
Outputs from all six controller states augment the three 
plant outputs, giving a total of nine feedback variables used 
in the variable-gain design example. The sampling period 
used is 0.31 seconds to correspond with an existing real- 
time simulation. 

There are several possible gain scheduled 
parameters, p, that might be used for the variable feedback 
gain matrix (7), but since this controller is for high-a 
research, a gain schedule that is proportional to a has been 
selected. The function used is 


where a is the reference value (5°) of the lowest a-trim 

case. Values of p for the four models used in the design 
arc 0, 2, 6 and 11. Each of the models have equal 
weighting in the global cost (9). 

I inear Analysis 

Table II gives the s-pianc, closed-loop eigenvalues and 
damping ratios for the four models. All closed-loop 
eigenvalues are computed in the z-planc (discrete), but are 
shown in the s-plane for simplicity. Each model 
corresponds to a flight condition as described in the 
previous section. The six most negative eigenvalues under 
model 1 correspond to actuator dynamics and the n z filter. 

The closed-loop short period eigenvalues arc at -3.19 
+ j 1 .91, with a damping ratio of 0.86. The two 
eigenvalues nearest the origin represent the phugoid which 
has a reduced frequency and is more heavily damped than 


Jc co|resp°ndin ? open-loop phugoid. The phugoid is not 
directly controlled, and the fact that the closed-loop 
phugoid is changed from the open-loop phugoid is due to 
coupling through the feedback measurements. The closed- 
7°P' ™P e , , l od eigenvalues for models 2 through 4 are 
-0.96 ± jl.17, -0.79 ± j0.84, and -0.72 ± jl.04 
respectively. Comparing all four models, the smallest 
dampmg ratio (excluding the phugoid and actuators) is 
U.57 which is the short period mode for model 4 The 
&°l da I ? pi ^ ra , ti0 for thc corresponding case is 
„„ ^I ablc A, 1 , 1 otJ } cr damping ratios are greater than 
U.6J. Thc controller (not shown) is completely stable 

except for one eigenvalue at the origin, due to thc a or a 
integrator for each model. 

Figure 1 contains the structured singular-value^ f or a 
multiplicative error when breaking the control loop for each 
Pk™ m P ut (toP P lo 0 or the plant output Gower 
plot). The numbers on the curves refer to the model Thc 
curve for model 1 is always above 1.0 for the plant input, 
and only goes below 1 .0 at thc lower frequencies for the 
plant output; the minimum is at 0.01 rad/scc, which is thc 
closed-loop phugoid frequency. Model 2 results show 
minimum singular values slightly less than 0.9 (excluding 
the phugoid frequency which is very lightly damped) at 
both the plant input and the plant output. The lowest 
singular value at the plant input is slightly greater than 0.6 
tor model 4, and the lowest singular value at the plant 
output is approximately 0.48 for models 3 and 4. This last 
aata indicates that a simultaneous diagonal perturbation at 
the plant output, as large as 48 percent, can insure stability 
at the frequency where the singular value is minimum. 

The maximum singular value for the loop transfer is 

vqWc fH fi i? f re 2 ‘ I he top plot conuins the singular 
values for all four models when the loop is broken at the 

plant input, and the bottom plot shows the singular values 
when the loop is broken at the plant output. The lightly- 
damped, open-loop phugoid frequency for models 1 to 3 
are shown by the sharp peaks in the singular- value curves, 
the sharp peak in the curves for model 4 represent thc 
open-loop short period which has a damping ratio of 

■ , where ‘) s the open -loop phugoid for model 4 is 
highly damped at 0.61. The plots show that the crossover 
frequency at the plant input is less than 9 rad/sec for 
model 1 and less than 3 rad/sec for thc other 3 models. 

All models show good attenuation characteristics at the 
high frequencies. 

Prior to analyzing the variable-gain controller in a 
nonlinear simulation, single-point designs were developed 
to compare linear and nonlinear simulation results for 
maneuvers with large anglc-of-attack excursions. 
k ^ u nce . h ? s S ^ 0wn that nonlinear simulation results 
should be similar to the linear simulation results, and gross 
differences between the two usually imply discrepancies in 
code or other anomalies in thc nonlinear simulation 
Hgure 3 contains linear simulation plots for a point design 

at a = 60°. This single-point design is different than thc 
frozen vanable-gain design point described previously for 
model 4. The only nonlinearity is a limiter that is inserted 
at thc plant input to prevent unrealistic control commands. 


The top portion of figure 3 shows curves for a, q, 0, n 
and y m . A ramp signal ,y m , with a slope of 1 rad/sec and 

a saturation of 1.0 rad is commanded. Thc a curve 
reaches thc commanded value in approximately 3.2 sec, 
and then settles out to a steady-stale error of approximately 


0.1 rad. The reason for this error is due to the saturation 
of 8 S which is shown in the lower portion of figure 3 

along with 8 y . Thc remaining control, 8 y , is normally 

commanded to zero degrees in steady-state, but 8 levels 
out at a new steady-state value to help compensate for the 
saturation of 8 s . Thc two integrator errors for a and 8 

leaving each variable offset from the 
commanded value. There is approximately seven more 

degrees of operating range for 8 y , and a pilot could use 

dus additional control to get a closer to the desired value 
An equivalent nonlinear simulation for this single-point 

"4S If nonhnear simulations of a variable-gain 
design, will be shown m the following section. 

ttQNLI NFAR simulation 

The airplane model and key features relating to the 
nonlinear batch simulation has been described previously. 
i™t ^ n in . C ? r emulation has three controllers for the 
longitudinal, lateral and directional axes, respectively The 
longitudinal controller was replaced initially with a fixed- 
gain controller based on a single-point design, and then the 
variable-gain longitudinal control law described in this 
nKnJn 6 ! 3teral “, d d * rec tional controllers remained 
“ n cou P lin S signal, normally included 

in the longitudinal controller to compensate for lateral - 
direcuona! rotation rates, is not included in the simulations 
. scn • • In thc non l'near simulation, variable-gain 

n C S ‘" g ,s . uscd f ° r all feedback gains, whereas linear 
interpolation is used for the feedforward gains. 


Practical Asn^ 

, J!} re ? Practical aspects that should be considered are 
described in this section. The first is the need to use multi- 
mode controllers and to smoothly change modes during 
light. As discussed previously under Design Example 
two modes have been included in the variable-gain design 
with a transition between the modes. 8 

SCC ° nd cons ' dc ration is to incorporate trim 
schedules m the nonlinear simulation. Linear designs are 
based upon flap and thrust levels for a particular trim case. 
In simulation, flap and thrust schedules are incorporated as 

a function of a to make the simulation perform as close as 
possible to the actual airplane. The schedules are shown in 
ligure 4 as a function of time since this data is taken from 

one of the simulation runs. The five curves are a leading 
and •railing-edge left flaps (DNL, DFL), power-lever angle 
command (PLAC), and thrust (TT). Units of the top four 
curves are degrees, and the bottom curve is in thousands of 
pounds. The thrust trim schedule is incorporated to 
simulate the acuon of a pilot who would normally push the 
throttle forward for increased thrust as the airplane pitches 
up. Thc open-loop thrust schedule is based upon the trim 

Lhrust (which is related to a) from the linear design. 
Notice that the thrust level lags the command bv 
approximately 1 second. y 


A third practical consideration is a deadband 
nonlinearity in the 8 y control. Figure 5 shows the effective 

mining angle as a function of thrust vectoring deflection 
and nozzle area at an altitude or 15,000 feet. As an 

example, for a nozzle area of 320 square inches and a 8 

V 


4 


deflection of 5°, there is no effective turning angle. The 
normal deadbands at the high-alpha cases considered are 

approximately ±6° out of a total range of ±25°. The 
approach used to help nullify this deadband effect is 
described later in this section. The thrust vectoring 
deadband has an effect on the control derivatives that are 

used for the linear designs. For cases where 8 y is trimmed 
to zero deflection, the deadband must be accounted for in 
all perturbations from trim. The 8 y control derivatives 

used in the linear designs described previously have been 
adjusted to compensate for this deadband. 

Single-Point Results 

Figures 6A and 6B illustrate nonlinear simulation 

results for the single-point design at a = 60°, and should 
be compared with the linear time responses in figure 3. 
One main difference in Figure 6 is that all angular units are 
in degrees. The command is a ramp input that starts at 
0.2 sec and reaches the desired value at 1 .2 sec. In figure 

6A the airplane is trimmed at a = 35°, and then 
commanded to reach an a of 60°. The a curve reaches 

60° in approximately 3 seconds, has a slight overshoot, 
and then settles to the desired value. The bottom two 
curves in the figure (5 S j and 5 y j ) represent the left 
stabilator and left thrust vectoring deflection controls. 
Although 8 S j is approximately at the saturation level, there 

is enough power to maintain a at the desired value. The 
deadband adjustment can be seen in the 5 y j curve. The 

effective thrust vectoring deadband is estimated as a 
function of nozzle area. This deadband is then added to the 

8 y | linear command to get the total 8 y command. As seen 

in the figure, 8 y is operating near zero deflection angle 
since it is not needed for this maneuver. 

In contrast, Figure 6B shows the same curves for an a- 

command of 90°. These results look more like the linear 
simulation results in Figure 3. The alpha curve overshoots 
90° and then settles at a steady-state error of 5° to 7°. 5 s j 

is in hard saturation and 8 y j is seen to peak at 15° and 

settle out at approximately 12°. There is much more thrust 
vectoring control left, and a pilot could easily use this 

control to increase a to the desired value. 

Variable-Gain Results 

Figures 7A and 7B illustrate the longitudinal variable- 
gain control using the four models described earlier in this 
paper. The curve labeled p represents the gain-scheduling 
parameter and the curve labeled ERROR represents the 
error between the command generator output, y , and the 

regulated variable (q or a). y m initially represents a 
commanded q of 0 deg/sec and then changes to 60 dcg/sec 
at 0.2 seconds. As a increases, the mode transitions to an 
a-command system and y m represents 60° of a. At 

8 seconds into the run, y m commands a to decrease to 


35°, and at 16 seconds, y m reverts back to the initial 
command mode of 0 deg/sec for q. 

The a curve starts at 5° trim and reaches 55° in less 
than 3.5 seconds, then slowly reaches the initial command 
of 60°. During this same time period the pitch attitude, 0, 
reaches 90°, with the flight path angle (not shown) peaking 
at 30° before dropping off. q reaches a peak of 40 
deg/sec. At 8 seconds, a starts to decrease to the 

commanded values of 35°, reaching this value within 
3 seconds with a slight overshoot. The Final command 
shows a decreasing to the initial trim value of 5° 
approximately 3 seconds after the command is initiated. 
The q curve settles out to 0 deg/sec in steady-state. Both 

controls, 8 g j and 8 y j, remain within the linear range. 

This contrasts with the single-point design where the 
stabilator saturates. With the exception of the initial 

transient response, 8 . operates near the deadband. The p 

curve has the shape of a, since the gain-schedule 
parameter is directly related to angle of attack. The 
ERROR curve represents the difference between the 
commanded value and the regulated variable, and is a 
function of the mode. Initially, ERROR is the difference 

between y m and q. As a increases, ERROR becomes the 
difference between y m and a. Finally, during the last 

8 seconds of the simulation, ERROR is again the 
difference between y m and q. Notice that ERROR goes to 

zero during all three steady-state periods. 

The final set of plots, figures 8A and 8B, illustrate the 
ability of the variable-gain design to hold a while 
performing a lateral maneuver. This maneuver has 
similarities to a Hcrbsl maneuver 13. The longitudinal 
control shows y m ramping from 0 deg/sec to 60° in 

2 seconds. Seven seconds into the run, when a is 60° 

and 9 slightly greater than 90°, a full lateral stick (PCA) is 
applied. The lateral command is maintained for 3 seconds 
as the airplane rolls about the velocity vector, illustrated by 

changes in Euler angles 9, <f>, and vp. At 1 1 seconds, PCA 
is returned to zero. At this time, there is a slight decrease 
in a, which then returns to the commanded value. After 
the rolling motion stops, the airplane nose pulls through 
negative 0 and then heads to nose level. At the end of this 
maneuver, die airplane is heading in the opposite direction, 
holding a at 60° and with a negative flight path angle of - 
70°. 

One important feature of the control approach can be 
observed from the curves § y (thrust vectoring symmetric 

command) and 8 y j. In the simulation, the longitudinal 

and lateral-directional thrust vectoring commands are 
summed prior to driving the control surfaces, with the 
longitudinal signal having priority. Excluding the initial 
transient while a is increasing, 8 y remains close to the 

deadband limits because of the integrator on (he 
longitudinal thrust vectoring control command. During the 
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lateral maneuver, full thrust vectoring moments are applied 
by the lateral-directional controller to maximize roll rate 

because 5 y c is maintained near zero deflection angle. A 

second benefit of maintaining 8 v near zero deflection 

angle is to reduce heating on the vanes by keeping them at 
uie edge of the exhaust plume. 

SUMMARY 

An optimal variable-gain output feedback control 
approach has been described and applied to a highly 
nonlinear flight regime that includes stall and post stall. 
The variable-gain formulation takes advantage of linear 
control design techniques to develop a nonlinear control 
law. In contrast to the standard industry approach of 
scheduling feedback gains using either an interpolation 
technique or a least squares fit, the variable-gain 
methodology is a control design approach where a global 
cost is minimized as a function of selected design 
parameters. Each of the operating points considered in the 
design are guaranteed to be stable. 

A discrete control structure is used in the application 
The structure allows proportional feedback, integral 
feedback, and incorporates additional filtering by feeding 
back the control position commands that go to the plant. 
Additional dynamic compensation is also included to help 
regulate the short period mode. The controller used for 
this application is sixth order. The equations implemented 
in the nonlinear batch simulation are incremental and use 
total measurable quantities. Advantages of this formulation 
include automatic trim and subtraction of slowly varying 
biases in the measurements. 

The airplane model has characteristics of a high 
performance, high angle-of-attack airplane. Thrust 
vectoring controls have been added to help stability 
augmentation and maneuverability in the post-stall flight 
regime. Aerodynamic data is available for the angle-of- 

attack range from -10° to 90°. Analysis of the open-loop 
short period mode indicates extremely low damping in the 
post-stall flight regime, and a few cases where the airplane 
is unstable. 

The design approach is applied to four flight conditions 
(ct = 5 , 15°, 35°, and 60°). A pitch rate command mode 
is used at low a, an angle-of-attack command mode is 

used at high a, with a transition region between the two 
modes. The design includes two controls (stabilator and 
llirust vectoring) and three measurements (q, a, and n ). 

Closed-loop eigenvalue analysis of all four operating 
conditions (excluding actuator dynamics and the phugoid 
which is not directly controlled) show that the lowest 

damping ratio is 0.57 for the short period mode at a = 
60°. 


Two maneuvers are illustrated for the variable-gain 
scheduling example. These maneuvers illustrate the 
smooth transition between modes, pitch-up and pitch- 
down capability, the ability of the controller to regulate at 
the commanded values with zero steady-state error, and the 
ability of the controller to hold a high a in the presence of 
lateral maneuvers. An important control feature is the 
ability to hold the thrust vectoring control near a 
commanded value (generally zero deflection), allowing 
maximum thrust vectoring capability for the lateral- 
directional system and helping to minimize heating of the 
vanes. 


REFERENCES 

t. Halyo, N.; and Broussard, J. R.: Investigation, 
Development, and Application of Optimal Output 
Feedback Theory - Volume I: A Convergent 
Algorithm for the Stochastic Infinite-Time Discrete 
Optimal Output Feedback Problem. NASA CR- 
3828, August 1984. 


Halyo, IN.; Mocrdcr, D. D.; Broussard, J. R.; and 
Taylor, D. B.: A Variable-Gain Output Feedback 
Control Design Methodology. NASA Contract 
NAS1-17493. NASA CR-4226, March 1989 


3. Halyo, N.: A Variable-Gain Output Feedback 

Control Design Approach. Presented at the AIAA 
Guidance, Navigation, and Control Conference 
Boston, MA. August 14-16, 1989 


4. Moerder, D. D.; Halyo, N.; Broussard, J. R.; and 

Caglayan, A. K.: Application of Precomputed 

Control Laws in a Reconfigurable Aircraft Flight 
Control System. ICS TM-86-102, Information and 
Control Systems, Inc., Hampton, VA. 1986. 

5. Ostroff, A. J.; and Hueschen, R. M • 
Reconfigurable Multivariable Control Law For 
Commercial Airplane Using A Direct Digital Output 
Feedback Design. NASA TM-85759, February 


6. Ostroff, A. J.; and Hueschen, R. M.: Investigation 
of Control Law Reconfigurations to Accommodate a 
Control Element Failure on a Commercial Airplane. 
Presented at the 1984 American Control Conference 
San Diego, CA. June 6-8, 1984. 

7 . Ostroff, A. J.: Techniques for Accommodating 
Control Effector Failures on a Mildly Statically 
Unstable Airplane. Presented at the 1985 American 
Control Conference, Boston, MA. 1985. 

8 . Broussard, J. R.: Design, Implementation and Flight 
Testing of PIF Autopilots for General Aviation 
Aircraft. NASA Contract NAS1-16303, Information 
and Control Systems Inc. NASA CR-3709, 1983. 


Comparison of linear and nonlinear time simulations 
for a single-point design at a=60° show very similar 
results. Similarities include rise-time, overshoot and offset 
errors when one of the controls is saturated. Experience 
has shown that nonlinear simulation results should be 
similar to the linear simulation results, and gross 
differences between the two usually imply discrepancies in 
code or other anomalies in the nonlinear simulation. 


9. Maybeck, P. S.: Stochastic Models, Estimation, and 
Control. Vol. 3, Academic Press, New York. 1982. 

10. Broussard, J. R.; and O'Brien, M. J.: Feedforward 
Control to Track the Output of a Forced Model 
IEEE Trans. Automatic Control, Vol. 25. August 
1980, pp. 851-854. 


6 



11. Broussard, J. R.; and Halyo, N.: Active Flutter 
Suppression Using Optimal Output Feedback Digital 
Controllers. NASA Contract NAS 1-16772, 
Information and Control Systems, Inc. NASA CR- 
165939, May 1982. 

12. Doyle, John: Analysis of Feedback Systems With 
Structured Uncertainties. IEEE Proc., Vol. 129, Pt. 
D, No. 6. November 1982, pp. 242-250. 

13. Herbst, W. B.: Future Fighter Technologies. 

Journal of Aircraft, Vol. 17. August 1980, pp 561- 
566. 


Table I - Short Peri od Open-Loon 


a 

Mach 

Eigenvalue 




Real 

Imag 

5 

.46 

-.55 

1.30 

10 

.33 

-.37 

.66 

15 

.28 

-.28 

.87 

20 

.26 

-.25 

.64 

25 

.24 

-.25 

.36 

30 

.22 

-.24 

.58 

35 

.20 

-.27 

.56 

40 

.19 

-.27 

.69 

45 

.19 

-.033 

.73 

50 

.19 

.019 

.83 

55 

.19 

.023 

1.17 

60 

.20 

-.017 

1.07 


Table II - 

Closed-Loon Eigenvalues 

Real 

MQBEL1 

Imag 

Damping 

-.622E+02 

.857E+02 

0.59 ■ 

-.622E+02 

-.857E+02 

0.59 

-496E+02 

.000E+00 

1.00 

-301E+02 

.OOOE+OO 

1.00 

-.166E+02 

.337E+02 

0.44 

-.166E+02 

-.337E+02 

0.44 

-.158E+02 

.000E+00 

LOO 

-.531E+01 

.600E+01 

0.66 

-.531E+01 

-.600E+0 1 

0.66 

-319E+01 

.191E+01 

0.86* 

-.319E+01 

—.191 E-+-0 1 

0.86* 

147E+01 

,47 8 E +00 

0.95 

-.147E+01 

-.478E+00 

0.95 

-.155E+00 

.OOOE+'OO 

1.00 

-.100E-01 

.739E-02 

0.80 

-.100E-01 

-.739E-02 

0.80 


MQPLL.2 


Real 

lmag 

-.621E+02 

.852E+02 

-.621E+02 

-.852E+02 

-.51 1E+02 

.000E+00 

-.306E+02 

.000E+00 

-.223E+02 

.000E+00 

-.153E+02 

.333E+02 

-.153E+02 

-.333E+02 

-.693E+01 

.000E+00 

-.461E+01 

.372E+01 

-.461E+01 

-.372E+01 

-.960E+00 

. 1 17E+01 

-.960E+00 

-.1 17E+01 

-.862E+00 

.000E+00 

-.187E+00 

.000E+00 

-.1 17E-01 

.133E+00 

-.1 17E-01 

-.133E+00 


MODEL 3 

Real 

Imag 

-.621E+02 

.851E+02 

-.621E+02 

-.851E+02 

-.517E+02 

.000E+00 

-.317E+02 

.OOOE+OO 

-.280E+02 

.000E+00 

-.150E+02 

.332E+02 

-.150E+02 

-.332E+02 

-498E+01 

.370E+01 

-.498E+01 

-.370E+01 

-.379E+01 

.675E+00 

-379E+01 

-.675E+00 

-.791E+00 

.841E+00 

-.791E+00 

-.841E+00 

-.929E-01 

.OOOE+OO 

-.566E-01 

.177E+00 

-566E-01 

-.177E+00 


MODEL 4 

Real 

Imag 

-.621E+02 

.849E+02 

-.621E+02 

-.849E+02 

-.507E+02 

.OOOE+OO 

-.329E+02 

.OOOE+OO 

-.285E+02 

.OOOE+OO 

-. 148E+02 

.333E+02 

-. 148E+02 

-.333E+02 

-.886E+01 

.106E+01 

-.886E+01 

-.106E+01 

-.221 E+01 

. 195E+01 

-.221E+01 

-. 195E+01 

-.720E+00 

.104 E+01 

-.720E+00 

-.104 E+01 

-.317E+00 

.OOOE+OO 

-.204E+00 

.432E-01 

-.204E+00 

-.432E-01 


* Short period 


Damping 

0.59 

0.59 

1.00 

1.00 

1.00 

0.42 

0.42 

1.00 

0.78 

0.78 

0.63* 

0,63* 

1.00 

1.00 

0.09 

0.09 


Damping 

0.59 

0.59 

1.00 

1.00 

1.00 

0.41 

0.41 

0.80 

0.80 

0.98 

0.98 

0.69* 

0.69* 

1.00 

0.30 

0.30 


Damping 

0.59 

0.59 

1.00 

1.00 

1.00 

0.41 

0.41 

0.99 

0.99 

0.75 

0.75 

0.57* 

0.57* 

1.00 

0.98 

0.98 
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Figure 1 Singular Value-Multiplicative Error. Figure 4.- Flap and Thrust Schedule. 
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Figure 2.* Singular Value-Loop Transfer. 
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Figure 5.- Effective Turning Angle as a Function of Thrust 
Vectoring Deflection and Nozzle Area. 



Figure 3.- Time Response for Design at 60°. Figure 6a - Two Independent Controls, Two Integrators, 

Alpha Command = 60. 
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